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Abstract. The neoclassical transport theory is applied to calculate electron cyclotron 
current drive (ECCD) efficiency in an axisymmetric tokamak in the low-collisionality 
regime. The tokamak ordering is used to obtain a system of equations that describe the 
dynamics of the plasma where the nonlinear ponderomotive (PM) force due to high- 
power RF waves is included. The PM force is produced around an electron cyclotron 
resonant surface at a specific poloidal location. The ECCD efficiency is analyzed in the 
cases of first and second harmonics (for different impinging angles of the RF waves) and 
it is validated using experimental parameter values from TCV and T-10 tokamaks. The 
results are in agreement with those obtained by means of Green's function techniques. 
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1. Introduction 

Electron cyclotron waves can efficiently drive a localized non-inductive current in 
toroidal devices for a number of applications. The main of these are found in the 
neoclassical tearing mode control pp, the fully non-inductive current drive in tokamaks 
[21 IS] and the bootstrap current compensation in stellerators 4\. ECCD results from the 
sensitive heating of electrons travelling in one direction in order to decrease their collision 
frequency, and thus enhance their contribution to the toroidal current, compared to their 
unheated counterparts moving in the opposite direction For an off-axis current drive, 
this current drive mechanism is offset by the mirror trapping of electrons in toroidal 
geometries that drives current in the reverse direction j^j. The ECCD efficiency is usually 
calculated through a bounce-averaged quasilinear Fokker-Planck treatment |H] . 

Electron cyclotron (EC) waves have recently attracted a great interest. Such 
waves exhibit the very important property of being able to be excited at a localized 
particular magnetic surface. This mechanism considers the introduction of EC waves 
at a minimum of the magnetic field, where the resonance condition v zr = (u — oose) /k z 
holds. Many tokamak experiments have reported that the ECCD efficiency decreases as 
the power deposition location is moved away from the plasma center, either by varying 
the magnetic field strength [§] or by changing the poloidal steering of the ECCD launcher 

ma- 

In the low power limit, the ECCD can be calculated from the relativistic, linearized 
Fokker-Planck equation using ray tracing codes jllj . If the effects of the radiofrequency 
(RF) quasilinear diffusion and the parallel electric field are included, the bounce- 
averaged, quasilinear Fokker-Planck codes can be used However, the nonlinearities 
associated to high power effects are not considered. In the present paper we analyze the 
ECCD efficiency in an axisymmetric tokamak, in the low collisionallity regime within 
the neoclassical transport theory. The ECCD is calculated including a ponderomotive 
force T . The tokamak ordering is used to obtain a system of equations that describe 
the dynamics of the plasma where the nonlinear ponderomotive (PM) force due to high- 
power RF waves is included. The PM force is produced around an electron cyclotron 
resonant surface at a specific poloidal location. The ECCD efficiency is analyzed in 
the cases of the first and second harmonics (for different impinging angles of the RF 
waves) and it is validated using experimental parameter values from TCV and T-10 
tokamaks. The results obtained are in agreement with those delivered by the linearized 
Fokker-Planck equation. 

2. Basic equations. 

Let us assume a plasma which contains only charged and neutral particles in a toroidal 
axysimmetric magnetic field. The hydrodynamic description of the plasma is taken in 
the neoclassical fluid approximation. In this approach, the continuity equation for the 
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averaging quantities respecting the RF field becomes 
On 

n + V • (?7 Q U Q ) = S na , 



dt 



where S na is the source term obeying the condition 

J2^s na = J^va [ ( -j^- ) rfv = o. 

a a J \ 01 / sQ 

where S na = f (df a /dt) s dv. For the moment equation, we have 

d(n a XJ a ) ^ 4=> „ ,__ tt tt \ , t-> (a) TJ (a) 



(1) 



(2) 



rn 



dt 



-V • P a - m Q V • (n„U a U a ) + + R^ j 



+ 9a 



n a E + -n a XJ a x B 

c 



(3) 



The subindex a refers to the particle species; m a , n a , and t/ a are the mass, electric 
charge, the density of particles and the velocity of the fluid, respectively. F a is the 
ponderomotive force, n a is the viscosity tensor and, finally, p a is the pressure defined 
by Pa — n aTai where T a is the plasma temperature. On the other hand, 

'df Q 



p(«) 
d(«) 



dt 



cO 



(it), 



(4) 
(5) 



sO 



is the force of friction between the particles of species a with neutrals. 

The system of equations and must be completed with the Maxwell equations 
for the averaging quantities 

V ■ E = Anq a n a0 , 

„ ^ 1<9B 
VxE = — — , 

c at 

VB =0, 



(6) 



_ _ IdE 4tt 

V x B = — — H g a n a0 U a . 

c at c 



By using an standard expansion with respect to the ratio of the gyroradius and the 
characteristic length, it follows that 

_ PLa 

£ ~ L ' 

where pL a is the Larmor radius and L is the characteristic length [T5j . 



3. Equations of zeroth and firsth orders. 

The system of equations CE}-© is reduced, to zero-th order terms, to 

V • n a0 U a0 = 0, 

V • P a = q a n a0 E + -n a0 U a0 x B 

L c 



(7) 
(8) 
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and the Maxwell equations become 



V-E = 


o. 


V x E = 


0. 


VB = 


o. 




47T 


V x B = 


— J«o- 




c 


ain that Eq 


= -V$ 



(9) 



In this case, the solution of the system ©-© has the form 

ut = o, 



+ 



9p a 



Ui=- 



d$ 

dip ' q a n a o dip 



1 dp Q 



R2 
n 



+ 



BqqB 



(10) 
(11) 

(12) 



where we have introduced the toroidal flux coordinates (ip, 8, £). Within this coordinate 
system, the contravariant forms of the magnetic field and of the fluid velocity for an 
axisymmetric tokamak are written as 

B = J(V)VC + VCx VV, 

U a = U d a e e + Ufa, 
where n a o = cte and the function A (ip) is unknown. 

Considering the inequality -Boc-^o ^ BqoBq and B\ 
equations (fTTj) and (fTi2j) to the form 

U 6 a = \{iP)Bl 



(13) 



-Boc-Bq we reduce the 



Ui 



JB b n 



d% 



+ 



1 dp a 



(14) 
(15) 



<9?/> g a n Q0 dip 

which are the zero-th order velocity equations. The corresponding toroidal current 
density is calculated from the relationship 



-JbJ 



+ 



1 <9p 



+ g«^ Q A (ip) Bq. 



(16) 



<9t/> g a ra Q0 dip 

Equations (0)-© containing terms to first order. Ignoring the source term they 
become 

dn aQ 



dt 



+ V • (n a0 V a ) = 0, 

= -Vp, 
+ q a 



d ( w tt0 u a ; 

dt 



V • 7f a - m a V • (n Q0 U a UJ + R^ 



n Q0 E + -^ a oU Q x B 

c 



(17) 



where P Q 
) , ■ 

L c0 



r»(a) I t»(«) 



Pa + 7Tq, Pa is the scalar pressure, Tt a is the viscosity tensor, R^*- n ) = 
s0 is the friction force from the particles of species a with neutrals and 
j — V • 7r a is the average ponderomotive force associated with the RF field 



acting on the particles. 
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4. Steady state equations. 

In a steady state, the above system of equations can be written in the form 

V • (VaoUa) = 0, 

m Q V ■ (VaoUaUa) = ~ Vp* ~ V ■ 7T a + 



+ q a 



n a0 E + -n a0 U a x B 

c 



(19) 
(20) 



where the ( component of the velocity becomes 



Ui 



c (d§ , 1 dp a \ B 0C B C Bi 



+ 
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JB b n 



(21) 



where /ig and are the poloidal and toroidal coefficients of viscosity, respectively. In 
the limit B gB® « Bq^Bq and B\ ps Bq^Bq, the toroidal current density is reduced to 



(0 



cn a q a 

HbJ 
+ 
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C/i C 



+ 
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(22) 



and 

A(V) 



(/i • B ) + m a r] aQ v an (B%) 



(B -F a )-k + 



2H_ 



d$ | 1 d Pa 
dip q a r] a0 dip 



+ m a i] a0 v an 



JB d 



d$ 1 dp 
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(23) 



5. The ponderomotive force. 

The RF ponderomotive force has several representations according to its functionality 
with respect to time. In this work, we chose the following expression of the time averaged 
ponderomotive force 



2 uj 



UJ 



UJ 



pa 



(24) 



where uj is the frequency of the RF wave, is the current density of particles of species 
a induced by the RF field E and uj^ a = Aim a q^/m a is the plasma frequency of particles 
of species a. 

Given the Ohm's law, we assume the conductivity tensor o~jk, which depends on 
the assumed characteristic frequency. The ponderomotive force is introduced thanks to 
a system of orthogonal coordinates (ei, e 2 , e 3 ) with the components of the conductivity 
tensor in the form 

UJ y/llV 



Cll = CT 2 2 



IUJ V 

4:71 1 — U 



Cl2 



"0-21 



AlTX 1 — 11 



(25) 
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Table 1. TCV tokamak data for the first and second harmonics. 





B 




a 


Ro 


P 


T 


<1 




(gauss) 


(cm -3 ) 


(cm) 


(cm) 


(MW) 


(KeV) 




1th harmonic 


1.43 x 10 4 


1.75 x 10 13 


25.0 


88.0 


1.0 


3.5 


10 


2nd harmonic 


1.43 x 10 4 


1.75 x 10 13 


25.0 


88.0 


1.0 


3.5 


10 



Table 2. T-10 tokamak data for the first and second harmonics. 
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Ro 


P 
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(gauss) 


(cm~ 3 ) 


(cm) 


(cm) 


(MW) 


(KeV) 




1th harmonic 


2.78 x 10 4 


0.54 x 10 13 


38.7 


1.5 x 10 2 


0.75 


6.3 


9 


2nd harmonic 


2.47 x 10 4 


0.54 x 10 13 


38.7 


1.5 x 10 2 


0.45 


3.8 


9 



where v = u 2 e /uj 2 ,u = Q 2 /u 2 and Q = eB/m e c. Here, we have considered the case of 
an extraordinary wave (E\, E 2 , 0). 

Thus, the corresponding components of the ponderomotive force take the form 

^_ 1 v Iraki / \ i 1 2 / n 

^i = 7r- t: \(l+u)\E\\ (26) 

07T (1 — U) 

^2 = —- 2 -t(l + u)\E\ 2 . (27) 

ovr (1 - u) 



(28) 



Now, by calculating the average value (B ■ T a ) assuming that 
B = B 9 e e + £ % 

J~ <x J~ ct1p^-1p~\~E a QCQ , 

and the equation, 

Bo ■ F a = FaoB®, (29) 
and finally substituting Po^ - (j2^j) in (}2T?j) and averaging, we obtain 



( '\ + u) \E\ 2 + Ay/ulm (E X E* 2 



\2 
U) 



B Imk2 
(Bo ■ T a ) = 2 

where it has been considered that 

d 2 uv 2 =d 2 ^ v{1 \ u) =0, 

\l-uf (l-uf 
and the Hamada coordinates [19J have been used. 

Finally, neglecting the attenuation of the RF wave, we obtain 

~- * Bnlmko v (1 + u) , ,o . . 

' b ^-'= a;, (1 -,/ w' < 3o) 

while the current density related to the ponderomotive force, from (|22jl . assuming the 
steady state, becomes 

■C(p) _ / T \ ~ n a g a /c 2 'l £ '| 2 + , 3 . 
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6. Analysis of results. 

In order to examine the expression for the current density (f3*T]) associated with the 
ponderomotive force, we adopt the criteria that the deposited energy by the RF wave 
has to be bigger that the internal enegy (E 2 > NT) so to include the nonlinear effects. 
This condition is satisfied on Tokamaks TCV ad T-10, where an analysis reported in 
PS HB] shows that E 2 ~ 5.75689 x 10 13 > NT = 98122.5 and E 2 ~ 3.37737 x 10 13 > 
NT = 54500.5 in these tokamaks, respectively. Such results indicate that the effect of 
the nonlinear ponderomotive force is highly important in determining the energy density 
introduced by the RF wave. We will consider a flux cylinder with a radius equal to the 
Larmor radius so to calculate this energy. 

The corresponding data for the TCV and T-10 Tokamaks are summarized in Tabled 
and Table El JH3 respectively. In both reports, the first and second harmonics of 
the EC waves were used provided that the introduction of the RF wave took place at 
the high magnetic field (HF) side. 

The power associated to the amplitude of the electric field follows from 
o 2P 

\E\ = , (32) 

UJEq 

where P is the power of the wave per time unit and Eq is the permittivity in vacuum, 
assumed to be a constant. 

For the imaginary part of the permittivity, 3/c2 = k 2 , we use the expression reported 
in J7j. In the case of the first harmonic, one has that 



where 



fe 2 ' 2y/n 



F(z) =- 
y ' 4 



P 2 t (2 - q) z" 2 



e 



qz 2 (5/2 - qf 1 



+ 



14 q \F 



2 



- + z + y/nzy/^ze z ($ - 1 



(33) 



(34) 



Here, $ (x) is the error function, flx = vtJc is the ratio of the thermal velocity to the 
light velocity in vacuum, and z = 2 (Q — uo) /0/3f.. 
Analogously, for the second harmonic, we have 

f = 51^,(1 + row-, (35, 



where 



_2*-l-,(l-l/2) 



2 2 - 1 - q 

The viscosity is neglected while the collision frequency between electrons and 
neutrals was taken in the form 

Ven = r- 1 = (3.5 x 10 4 T 3 / 2 /n e ) 1 . (36) 

It is important to notice that the current density is highly unstable and it depends 
strongly on the Q 2 /uj 2 relationship. However, it is possible to find an interval where 
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the current density stabilizes and, furthermore, its values reproduce those experimental 
ones reported in [THJ EH] , as can be observed in figure [T] and figure El 

From figure El and figure 0J we observe a process in which the Fish and Ohkawa 
mechanisms weaken each other, as reported in (7j. The general behavior is in good 
agreement with that described in [HI CHI, considering that their calculation was obtained 
from the linearized Fokker-Planck equation. 

Here, the density profile has been modelled as a parabolic one in order to analyze 
the current, n e = n 0e ^1 — | O^jj , where a is the Tokamak minor radius. 

It can be noticed in figure El that, to first order terms in the parameter e, the driven 
current density increases with the radius. 

7. Conclusions. 

The development of a driven current density expression that takes into account the 
ponderomotive force created by EC waves, has required the use of the neoclassical 
transport equations up to first order terms with respect to the parameter e = p/L at a 
steady state. 

The driven current density has been initially obtained in a system of toroidal flux 
coordinates. That description of the current density is transformed in terms of the 
Hamada coordinates ^H] which is necessary for its validation with experimental results. 
Thus, the expression for the ponderomotive force reported in is written in a local 
system of orthogonal coordinates (ei, e2, £3), where e3 is parallel to the toroidal magnetic 
field. 

The driven current density generated by an extraordinary wave at the cyclotron 
resonance of electrons, is analyzed as a function of the Q 2 /u 2 ratio. This is accomplished 
by using the parameters of the TCV and T-10 Tokamaks at the first and second 
harmonics, assuming that the introduction of the wave takes place at the HF side. 
From this results, we have obtained an interval of frequencies, in agreement with the 
experiments, where the current shows a stable behavior. 

In the particular case of a parabolic profile, it has been shown that the ECCD 
increases with the radius, at first approximation. Finally, it is important to notice that, 
according to [Hj, the efficiency is higher for the second harmonic as it is shown in 
figure 
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Figure 1. Current efficiency plotted against £1 2 /oj 2 for the first harmonic with 
parametric values from the tokamaks a) T-10 with / = 140/1.06 Ghz (solid) and 
/ = 140/1.12 Ghz (dashed), and b) TCV with / = 82.7/1.484 Ghz (solid) and 
/ = 82.7/1.5 Ghz (dashed) 



Figure 2. Current efficiency plotted against e for the first harmonic with values taken 
by parameters of the tokamaks a) T-10 with r = 0.0 cm (solid) and r = 88 cm (dashed), 
and b) TCV with r = 0.0 cm (solid) and r = 88 cm (dashed). 



Figure 3. Current efficiency plotted against £l 2 /u> 2 for the second harmonic with 
parameters from the tokamaks a) T-10 with / = 140 x 0.7 Ghz (solid) and / = 
140 x 0.715 Ghz (dashed) and b) TCV with / = 82.7 x 1.331 Ghz (solid) and 
/ = 82.7 x 1.339 Ghz (dashed). 



Figure 4. Current efficiency plotted against e for the second harmonic with parameter 
values from the tokamaks a) T-10 with r = 0.0 cm (solid) and r — 88 cm (dashed), 
and b) TCV with r = 0.0 cm (solid) and r — 88 cm (dashed). 



Figure 5. Current efficiency plotted against the minor radius r for the second 
harmonic with parameters from the tokamaks a) T-10 and b) TCV. 
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